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ABSTRACT 

We have developed a model for gamma ray emission in jets of active galactic nuclei 
in which particle acceleration takes place at a shock in the relativistic jet plasma due 
to a massive star in the central region of the host galaxy moving through the jet. The 
gamma rays are produced in a pair-Compton cascade in the radiation field of the star 
initiated by accelerated electrons. Our model may account for the observed GeV to 
TeV gamma ray spectrum and variability of Markarian 421 and other blazars detected 
by the EGRET instrument on the Compton Gamma Ray Observatory. 
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1 INTRODUCTION 

In recent years many blazar type active galactic nuclei 
(AGN) have been detected in high energy 7-rays by the 
EGRET telescope (von Montigny et al. 1995, Thompson et 
al. 1995). Of these, The BL Lac object Markarian 421 has 
been observed at TeV energies by the Whipple Observatory 
(e.g. Punch et al. 1992) and HEGRA (Petry et al. 1996). 
A further BL Lac object not detected by EGRET, Markar- 
ian 501, has also been detected by the Whipple Observa- 
tory (Quinn et al. 1995). The 7-ray emission from blazars 
is strongly variable on different time scales (from less than 
an hour to weeks) in different energy ranges (see Gaidos et 
al. 1996, Mattox et al. 1997). In the case of the Markarian 
421, the X-ray variability is correlated with TeV variability 
(Takahashi et al. 1996a, 1996b; Buckley et al. 1996) with 
a time scale of ~ 1 week. Interestingly, the X-ray emission 
from this source seems to vary almost periodically on a time 
scale of ~1 day (Takahashi et al. 1996a). Such variability is 
not naturally expected in the present popular scenario for 
the production of 7-rays in a blob moving relativisticly along 
the jet in a blazar. 

Recently, Dar and Laor (1996) proposed that 7-ray pro- 
duction in AGN jets may be due to the collision of small 
clouds with radii ^ 10^^ cm and densities ~ 10^^ cm~^ with 
a highly coUimated relativistic proton beam with 7 ~ 10*. In 
the present paper we consider another scenario as a possible 
explanation of coUimated 7-ray production in blazars. 

In the central regions of active galaxies the processes of 
star formation are probably very efficient. For example, in 
M32 the central stellar density exceeds 10^ M0 pc~'^ (Lauer 
et al. 1992). It seems obvious that many stars must collide 
with the jet plasma. Here we investigate the consequences 
of such frequent stellar wind - jet plasma collisions for the 



production of highly coUimated 7-ray beams in blazars. The 
general scenario is illustrated by Fig. 1(a). 



2 STAR COLLIDING WITH JET PLASMA 

Let us consider a very simple scenario in which a pow- 
erful jet, propagating from the central engine of active 
galactic nuclei, meets a massive star. The jet may consist 
of relativistic electron-proton (neutral) plasma with power 



of injection of protons with Lorentz factor 7 ~ 10 into the 
jet, and nip, me are the proton and electron masses. The 
kinetic powers of jets in radio-loud AGNs may be up to at 
least ~ 10*^ erg s~^ (Rawlings & Saunders 1991, Celloti 
& Fabian 1993, Falcke & Biermann 1995). Assuming the jet 
propagates within a cone with opening angle 6, at a distance 
I along the jet the ram pressure of the jet plasma is 



erg cm , (1) 

where L46 is the jet power in units of 10*® erg s~^, 65 is 
its opening angle in units of 5°, and h is the distance in 
parsecs. 

Stars in the central regions of galaxies move with high 
velocities, of the order of (3 — 6) x 10"^ km s~^ , and are char- 
acterised by very intense winds and strong surface magnetic 
fields. For example the mass loss rate of Wolf-Rayet (WR) 
stars and the terminal velocity of the matter outflow are 
MwR ~ (0.8 - 8) X lO-^M© yr-\ and ~ (1 - 5) x 10=* 
km s~^ (Lang 1991). Young massive OB stars also have in- 
tense stellar winds with Mob ~ lO'^M© yr~*, and w2P ~ 
(1 - 3) X 10^ km s"^ (e.g. Lang 1991). 

Some of these massive stars will obviously cross the jet 
cone from time to time. When this happens, the interaction 
of relativistic plasma in the jet with the stellar wind (and the 
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Figure 1. Schematic representation of the interaction of a star with an AGN jet. (a) The geometry (not to scale) showing stars orbiting 
around the centre of the galaxy and showing a star entering the jet, having opening angle 8, a distance I above the central engine, (b) 
Expanded view of the region around the star showing the stellar wind and the relativistic jet plasma interacting to produce a double 
shock and contact discontinuity (CD). In region A the jet plasma shock is quasi-parallel, and in region B is quasi-perpendicular. Particles 
are accelerated efficiently at the relativistic shock in the jet plasma, and in region B may be highly coUimated along the magnetic field 
direction which is expected to be parallel to the direction of jet plasma flow (see text). 



associated stellar magnetosphere) will result in the forma- 
tion of shock waves similar to those discussed in the context 
of collisions of winds in binary systems involving early type 
stars (e.g. Eichler & Usov 1993). In the present context, the 
location of the shocks is determined by the parameters of 
the jet and the star. A double shock structure, with a con- 
tact discontinuity between them, should form (see Fig. lb). 
The shock in the stellar wind will be non-relativistic while 
that in the jet plasma will be relativistic. 

In principle the jet pressure can be balanced by the 
stellar wind pressure or the stellar magnetic field pressure. 
For stars with a dipole magnetic field and a strong wind, 
the structure and strength of the magnetic field around the 
star is (Weber & Davis 1967, Usov & Melrose 1992) 

( {R/rf, for i? < r < (dipole), 
B{r) ^BsX I C\R/rf, for < r < Rr]-^ (radial), (2) 
I ■qC^{R/r), for r > Rr)'^ (toroidal). 

where ^ = rA/R, ta is the Alfven radius which depends on 
the star's parameters and is in the range ~(1-3)R, where R 
is the radius of the star (Usov & Melrose 1992), and rj = 
"Urot/vcxj where Vrot is the surface rotational velocity of the 
star. Typically for early-type stars v^ot = (0.1 — 0.2)iioo (e.g. 
Conti & Ebbets 1977). 

If the stellar wind is very strong, the jet pressure is 
balanced by the wind pressure, 

Pwind = ~ 1.6 X 10*M_5«3iir2'(r/i?)"'erg cm-^ (3) 

where the mass loss rate of the star is M = M-^Mq 
yr~^, the wind velocity is Vac = 3 x lO^vs cm s~^, and 
the radius of a star is R = 10^'^-Ri2 cm. The double shock 
structure is then located at a distance rsh from the star given 
by 

r,h/i? ^ 10iM]^^vy^R^^eshL-J^\ (4) 
The maximum power which may in principle be extracted 



from the jet by the shock can be estimated from 

L,h = Ljir,i,/9lf ^ 1.5 X 10^^M_5^;3 erg s'\ (5) 

Note that this depends only on the stellar wind power, but 
not on the jet power or distance from the central engine. 

If the jet pressure can not be balanced by the stellar 
wind pressure, it may be balanced by the pressure associated 
with the star's mag netic field, B{rf/S-K, provided 

B'i > 0.6M^Svl^^R^2^ (6) 

where B3 is Bs in units of 10^ G. Then the shock location 
would be given by 

r^h/i? ~ 5.5(B3e5/ii46''")''''- (7) 

In this case, the maximum power that can be extracted by 
the shock is 

Leh ~ 4.4 X 10''i??2(B3L46e5-'zr')'^'- (8) 

In the most extreme case of the jet pressure dominating 
over the wind and magnetic pressures of the star, the jet col- 
lides directly with the stellar atmosphere. The jet pressure 
dominates over the magnetic pressure if 

L46 > 2.7 X lO^Bielll, (9) 

and the jet pressure dominates over the wind pressure if 

L46 > WM-^vsOlllR^i. (10) 

If both conditions are met, the jet initiates a very strong 
outflow of matter from the star, which pressure balances 
the jet pressure very close to the star's surface. 



3 ACCELERATION OF PARTICLES 

Electrons and protons are expected to be accelerated by the 
first order Fermi acceleration mechanism (e.g. see reviews 
by Blandford & Eichler 1987, Jones & Ellison 1991). If the 
diffusion coefficient is proportional to energy we can write 
the acceleration rate as 
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E = xZecB erg (11) 

where B is in gauss, Ze is in statcoulombs, and x depends on 
the details of the acceleration mechanism. For acceleration 
at a shock with velocity 0.1c values of x as high as 0.04 
or 1.6 X lO""* are possible if the shock is perpendicular or 
parallel respectively (e.g. Protheroe 1997). 

For typical parameters of massive stars (surface tem- 
perature of the order of a few 10* K, and other parameters 
mentioned above, Lang 1991), accelerated protons are not 
expected to encounter sufficient target matter or radiation 
for interactions in the shock region. The interaction lengths 
(in cm) for protons in the wind plasma and stellar radiation 
are approximately 

\pp^3x lO'^MZlvzRl^irA./Rf, (12) 

Ap^ ^ 3.4 X 10'''T4-^(r,h/i^)^ (13) 

where the star's surface temperature is taken to be T = 
10*74 K. These interaction lengths are longer than the char- 
acteristic distance scale, i.e. the shock radius. We therefore 
concentrate on electron acceleration in this paper. Using the 
standard formula for the synchrotron energy loss rate one 
obtains, in the absence of other losses, the maximum elec- 
tron energy 

ET" = 6 X 10*x'^^-B-i/2 GeV. (14) 

As shown in Section 2, collisions of stars with difi'erent 
parameters (M, Vod, B, I) with jets having difi'erent parame- 
ters {Lj,6) will result in the formation of shocks at various 
distances from the star. Since the shock parameters will be 
different for the shock in the stellar wind and the shock in 
the jet plasma, which we shall refer to as "the wind shock" 
and "the jet shock" , we shall discuss them separately below. 

3.1 Shock in the jet plasma 

The jet shock may be very effective for particle acceleration 
because it occurs in the relativistic jet plasma. Observations 
of the inner parts of strong jets show that the magnetic field 
is mainly directed along the motion of the jet plasma (Saikia 
& Salter 1988). In such cases the region of the shock near- 
est to the central engine (region A to the left of the dashed 
line in Fig. 1) is quasi-parallel, while the region of the shock 
farthest from the central engine (region B in Fig. 1) is quasi- 
perpendicular. The strength of the longitudinal component 
of the magnetic field in the shock depends on the total mag- 
netic fiux, $, through the jet. At large distances from the 
base of the jet [61 ^ rf , where rf is the jet radius at the 
base of the jet), the magnetic field scales inversely with the 
jet cross-sectional area 

Bj ^ 1.5 X 10-^$326'5"'«r', (15) 

where the $32 is the magnetic flux in units of 10^^ G cm^. 

The section of the shock farthest from the central engine 
(region B in Fig. 1) seems to be the most promising location 
for the production of directly observable fluxes of 7-rays in 
the scenario discussed here. The reason being that in this 
region the relativistic shock is quasi-perpendicular and is 
likely to accelerate particles which are strongly coUimated in 
the shock surface (see e.g.. Kirk & Heavens 1989, Ostrowski 
1991). Particles can be also accelerated almost rectilinearly 



by shock drift acceleration (Jokipii 1987) discussed by Begel- 
man & Kirk (1990). The degree of particle coUimation can 
be of the order of 1/7||, where 7|| is the Lorentz factor for 
the velocity of the shock measured along the direction of the 
magnetic field lines. A very strong anisotropy of accelerated 
particles is then possible provided that the magnetic field is 
highly ordered, with no strong turbulence present. 

The distance of the shock from the star is given by 
Eq. (^) and the magnetic field at the shock can be esti- 
mated from Eq. (|l^). This enables us to estimate the ac- 
celeration rate and the maximum electron energy allowed 
by synchrotron losses. Using Eq. (^^, and Eq. ( pq ) for the 
magnetic field in the jet, one finds 

ET" = 1.5 X 10*'x'''^*M '"'^5/1 GeV. (16) 

For example, for x > 0.04 (expected for a perpendicular rel- 
ativistic shock) and h = 0.05 one obtains E'^'^'^ > 1.6 x 10* 
GeV. Higher energies are possible for a highly anisotropic 
particle distribution because energy losses by synchrotron 
radiation depend on [B sin a)^ and sin a can be small. 

3.2 Shock in the stellar wind 

In this case the shock may form as a consequence of a bal- 
ance between the jet pressure and the stellar wind pressure, 
or the magnetic pressure, or as a result of a direct colli- 
sion of the jet plasma with the stellar atmosphere. The lo- 
cation of the shock in first case is determined by Eq. (^. 
If ''sh ^ Rv~^ J then the shock is quasi-perpendicular be- 
cause it is formed in the toroidal magnetic field region. If 
^Rri~^ then the shock should be quasi-parallel since it 
is formed in the radial magnetic field region. 

Unless the wind shock is formed very far from the 
star the shock will be non- relativistic {voc ~ 0.01c), quasi- 
parallel, and the magnetic field in the stellar wind will be 
higher than in the jet plasma. Hence, the maximum energy 
will be much lower than for the jet shock, and is likely to 
be determined by synchrotron losses rather than by inverse 
Compton scattering. It is therefore likely that electrons ac- 
celerated at the wind shock may be responsible for the X-ray 
emission as a result of synchrotron radiation. To see whether 
this is reasonable, we consider the energy of synchrotron 
photons emitted by the highest energy electrons. 

Ex ~ rrieC (B/B (17) 

where B^ = 4.4 x lO" G, and 7max = E^"""^ /nieC^ (given 
by Eq. UM. For Ex ~ 10 keV one obtains an acceleration 
rate parameter x ~ 10~^ which is reasonable for a parallel 
non-relativistic shock. 

It is probable that the accelerated electrons may follow 
magnetic field lines to region B which are likely to become 
directed along the jet axis. As the electrons in this region 
cool, their pitch angles will be decrease such that their X-ray 
emission will be preferentially directed in the jet direction. 



4 PAIR-COMPTON CASCADE 

We consider the pair-Compton cascade initiated by relativis- 
tic electrons accelerated at the jet shock, such that they 
travel along the shock front, to estimate the spectrum of 7- 
rays emerging from the AGN in the present scenario. Since 
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Figure 2. Shock radius divided by mean free path of electrons for 
ICS (soUd curves), and shock radius divided by mean free path 
for photon-photon pair production (dotted curves), at r = r^ij 
multipled by {r^^^/R). Lower curves correspond to T = 3 X 10* K 
and R = 10^^ cm; upper curves correspond to T = 10^ K and 
R = 2x 10" cm. 



the shock is quasi-perpendicular and relativistic we consider 
injection of electrons with a power law spectrum index less 
than or equal to 2. The cut-off in the electron spectrum will 
be determined by the balance between particle energy gains 
from shock acceleration and losses by inverse Compton scat- 
tering, and can reach 10* GeV or higher as discussed below. 

We shall consider cascades in the radiation fields of OB 
stars for which we adopt T — 3x 10* K and R = 10^^ cm, and 
cascades in the radiation fields of Wolf-Rayet stars for which 
we adopt T = 10^ K and E = 2 x lO" cm. We have devel- 
oped a pair-Compton cascade code to calculate the emerging 
gamma-ray spectrum for the case of injection of electrons at 
the shock radius at the boundary between regions A and B 
and travelling in the jet direction (see Fig. lb). This code 
takes account of the anisotropic nature of the radiation field 
and is based on earlier work (Protheroe 1986, Protheroe et 
al. 1992, Protheroe and Biermann 1996). Using this code, 
we have calculated the ratio of the shock radius (character- 
istic distance scale for this problem) to the mean free path 
of electrons for inverse Compton scattering, and the ratio of 
the shock radius to the mean free path for photon-photon 
pair production in the anisotropic radiation field at the in- 
jection point. This ratio, multiplied by (vsh/R), is plotted in 
Fig. 2 and we find that the shock radius equals the mean free 
path for electrons of energy 10* GeV at the injection point 
if rsh = 5.2R for OB stars, or Vsh = 147? for Wolf-Rayet 
stars. Using Eq.(^) we could obtain these shock radii with, 
for example, the following parameters: M_5 = 0.1, — 1, 
L46 = 0.1, and k = 0.05 (OB stars); Ms = 8, us = 1, 
L46 = 30, and h = 0.05 (WR stars). Clearly, these are per- 
fectly reasonable sets of parameters which would give rise to 
pair-Compton cascading by energetic electron injection. 

Finally we check whether 10* GeV is a reasonable max- 
imum energy based on the acceleration rate and the energy 
loss distance for inverse Compton scattering (given approxi- 
mately by the mean free path in the Klein-Nishina regime) . 
Using Eqs. (^ and (^ the acceleration distance is given 

by 



Kacc = 2.2 X 10\-^E,'i>s2^elll cm, (18) 

where Ee is in GeV, giving x^cc = 1.4 x 10^^ cm for x = 0.04, 
'1'32 = 1, ^5 = 1 and h = 0.05. This acceleration distance 
is comparable to, but somewhat lower than, the shock radii 
estimated above confirming that 10* GeV is a reasonable 
maximum energy and that the maximum energy is probably 
determined by inverse Compton scattering. 

We have calculated the emerging gamma-ray spectrum 
for a power-law injection spectrum of the form 

( V n no^ 

IK = llGiV j ^'^^ (^'^^ 

for the standard shock acceleration spectrum, a = 2, and 
for a = 1.5 representing the fiatter spectra that are possible 
for acceleration at relativistic shocks (Ellison et al. 1990). 
Results are shown in Fig. 3(a) for OB stars and in Fig. 3(b) 
for Wolf-Rayet stars for a range of shock radii spanning those 
discussed above. Specifically we show results for Vsh/R = 1.7 
and Tsh/R = 17 for OB stars, and fsi-i/ R — 4.7 and rs\i/ R = 
47 for Wolf-Rayet stars. 

From Fig. 3 we see that when the shock is located far 
from the star, i.e. a;J^t(10* GeV) ^rgh, i.e., for less powerful 
jets and/or large I, significant gamma-ray emission will occur 
at TeV energies. Conversely, if the shock is near the star, 
i.e. x?2;(10* GeV) <C )"ah, the 7-ray spectrum is significantly 
steepened by photon-photon pair production in the stellar 
radiation field. For a fiat electron spectrum the TeV 7-ray 
fiux can be comparable to or dominate the GeV 7-ray fiux 
(see upper dashed lines in Fig. 3) as observed during TeV 
7-ray outbursts from Mrk 421 (Buckley et al. 1996). 



5 DISCUSSION 

The expected 7-ray luminosities (if recalculated for the 
isotropic case) are 

L.^ « 4 X lO^^iLshCfla, (20) 

where Lsh is given by Eq. (^) or (|^, /i is the efficiency of 
conversion of jet power into gamma rays at the shock, and 
a — 10~^Qf_3 rad is the angle of the cone of 7-ray emission. 
We assume that the axis of the emission cone is directed 
outwards from the source of the jet, and so the cone sweeps 
across the sky as the star crosses the jet. The expected vari- 
ability time scale associated with the line of sight to the 
observer passing through the cone of emission is then 

tvar ~ W asVi:^h, (21) 

where = 3 x loci's, 3 cm s~^ is the transverse velocity of the 
star. For a^z ~ l,'i's,3 = 1, and h = 0.05 this gives tvar ~ 6 
days, which is comparable with the activity period observed 
simultaneously from Mrk 421 in X-rays (Takahashi et al. 
1996a) and TeV 7-rays (Buckley et al. 1996) . It is interesting 
that Eqs. (|2^) and ( pl| ) predict that for a narrower emission 
cone 7-ray fluxes should be higher and variability time scales 
should be shorter. 

A distant observer may see emission modulated with the 
rotational period of the star if the star's magnetic axis does 
not coincide with its rotational axis. In this case, the wind 
and magnetic field strengths and structures should vary al- 
most periodically in the region of the shock with time scale 
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log(£'/GeV) log(£'/GeV) 

Figure 3. Emerging 7-ray spectra resulting from injection of a power-law electron spectra of the form E~^'^ (upper curves) and 
(lower curves) extending up to 10* GeV at the shock radius in a direction tangential to the shock, (a) T = 3 X 10* K, i? = 10^^, 
^ah/R = 1-7 (solid curves) and r^Y^/R = 17 (dashed curves); (b) T = IQS K, R = 2 X 10", r^hZ-R = 4.7 (solid curves) and r^h/R = 47 
(dashed curves). 



given by 



trot = 2ttRv-1 ^ (1 - 2) X i?i2«3"' days, (22) 

where Vrot = (0.1 - 0.2)«oo = 3 x 10*(0.1 - 0.2)v3 cm s"^ 
is the star's rotational velocity. This period is comparable 
to the possible quasi-periodic variability detected from Mrk 
421 (Takahashi et al. 1996a). Note also that winds from early 
type stars are very unstable, and the jet plasma may contain 
some irregularities. This may result in a slight change in the 
location of the shock which may significantly change the 
direction of a narrow cone in which 7-rays are emitted, and 
may be the reason for the flickering of the 7-ray emission on 
a time scale of fractions of an hour (Gaidos et al. 1996). 

Many stars are likely to be emerged in the jet at any 
one time. For a stellar density of ~ 2 x 10^ M0 pc~'^ within 
~ 1 pc (observed in M32, Lauer et al. 1992), about ~ 40 
stars with an average mass of ~ 10 M0 would be found 
inside a jet with opening angle 5° within ~ 1 pc of the 
central engine. The presence of many stars inside the jet 
significantly increases the probability of a distant observer 
being inside the 7-ray emission cone of one of the stars. 

In conclusion, we have developed a model for gamma 
ray emission in AGN jets in which particle acceleration takes 
place at a shock in the relativistic jet plasma produced as 
a result of a massive star in the central region of the host 
galaxy moving through the jet. The 7-rays are produced in 
a pair-Compton cascade initiated by accelerated electrons in 
the radiation field of the star. Our model may account for 
the observed GeV to TeV 7-ray spectrum and variability of 
Markarian 421 and other EGRET blazars. 
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